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Приближенное вьїчисление двойньх интегралов 
от бььстроосциллирующих функций с использованиєм 
лагранжевой полиномиальной интерлинации 


У статті досліджуються кубатурні формули обчислення подвійних інтегралів від швидкоосцилюючих 
функцій з використанням лагранжевої поліноміальної інтерлінації на класі диференційовних функцій. 
Інформація про функцію задана її слідами на системі взаємно перпендикулярних прямих. Отримані 
оцінки похибки наближення. 

Ключові слова: інтеграли від швидкоосцилюючих функцій двох змінних, 

кубатурні формули, лагранжева поліноміальна інтерлінація функцій. 


Сифбагиге Гогтиїіа58 ої арргохітаїе саїсшіацоп ої адоцйбіе іпіезгаїя ої рієрІу озсШаюогу Гипспоп5мій ицзе фе 
І аєтапее роїупотіа! іпіегійпеайоп ої Гипсйопя аї Ше сіаз5 ої дфіНегепцабіе Гипспопя аге іпуезіївакед їп Фе 
агісіег. Паїогтайоп абойі Гипсйоп 15 5ес ої уаїше58 оп (ре орійта! регрепаїсиіаг Ппез5. ТБе е5(ітайопя ої 
егтог ої арргоасріпє їо Ше сибакиге їогтиіаз аге ргезепіеа. 

Кеу УУ/огаз: іпіеотаї8 ої рієрІу о5сіШагогу Гапсйоп5, сибагиге Гогтпиіаз, І аєгапєе роїупопааї 
іпіегГіпеацоп ої Гипспоп5. 


В статье исследуются кубатурнье формульт приближенного вьтчисления интегралов от бьтстроосцилли- 
рующих функций двух переменньх с использованием лагранжевой полиномиальной интерлинации на 
классе дифференцируемьїжх функций. Информация о функции задана еб следами на системе оптимально 
вьбранньх взаймно перпендикулярньх линий. Полученьт оценки погрешности кубатурньх формул. 
Ключевьге слова: интегральт от бьюстроосциллирующих функций двух переменньтх, 
кубатурньєе формульі, лагранжевая полиномиальная интерлинация функций. 


Вступ 


Задача наближеного обчислення подвійних інтегралів від швидкоосцилюючих 
функцій є однією з найбільш важливих задач цифрової обробки сигналів. На даний 
час виникає необхідність наближеного обчислення таких інтегралів за допомогою 
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інформаційних операторів різних типів. Як дані можуть бути значення функції у вузлових 
точках, сліди функцій на системі взаємно перпендикулярних ліній тощо. Таку задачу 
ефективно дозволяє розв'язувати апарат інтерлінації функцій | 1), (21. В |3) вперше роз- 
глядалась задача наближеного обчислення інтегралів від шидкоосцилюючих функцій у 
випадку, коли інформація про функцію задана значеннями у вузлових точках, а в |4|, |3| 
представлена задача наближеного обчислення 2П1)-коефіцієнтів Фур'є за допомогою 
інтерлінації функцій у випадку, коли інформація про /(х, у) задається слідами на 
системі взаємно перпендикулярних прямих. Більш детально з такими кубатурними 
формулами, а також з їх тестуванням, можна познайомитись в (6), (7. 

Мета даної статті- дослідити питання оптимального вибору вузлів та ліній для 
запропонованих кубатурних формул. 

В |З) вперше було згадано про обчислення подвійних інтегралів від швидко 
осцилюючих функцій як повторних, з оптимальним вибором, вузлів. Питання побудови 
кубатурних формул для наближення подвійних інтегралів від швидкоосцилюючих 
функцій у випадку, коли інформація про функцію /(х, у), задана слідами на системі 


оптимально вибраних взаємно перпендикулярних прямих, розглядається в даній статті 
вперше. 
Постановка задачі: побудувати кубатурну формулу наближеного обчислення 

інтеграла 

11 

І, о)- Ї Г.С мо )зіпох, зіпохахах, 

-1-1 
з використанням лагранжевої поліноміальної інтерлінації функцій на класі дійсних 
функцій, визначених на ва|-мії, і таких, що | Ї (Ризра) (ж, у|«м . Отримати оцінку 


похибки кубатурної формули. 
1 Найкраща в і бів 4д 2 2,1,2 лагранжева поліноміальна інтерлінація. 
Наведемо деякі теореми. 
Теорема 1. |2| Нехай 2(г) є С"(І), І «0,11, Іх г «п, І, ,є(1)- інтерполяційний по- 
ліном Лагранжа степеня п- 1 функції 2(г), 
С Р 


І 187) то Ук дини (1), ва (т) - П 
Ка! 


ітінк р Й 


Кеїп 


із властивостями 
І леПдвє()кКаіп, -охахі чі З..ХІ,іХіХхрХо0, 
Тоді для залишку К, (г): 2(1)- І, ,2(г) справедливе інтегральне зображення 


п - Со 
ваФеУ |в Уа 
о с) 


При знаходженні найкращої в В; ЕЕ й 4д 2 ,1,2 лагранжевої поліноміальної 


інтерлінації на системі взаємно перпендикулярних прямих слід враховувати, що: 

1) залишок наближення інтерлінації дорівнює операторному добуткові залишків по 
кожній змінній окремо; 

2) поліноми степеня п з найменшим відхиленням від нуля в метриці ва - 


це поліноми з коефіцієнтом одиниця при старшому степені, що є розв'язанням екс- 
тремальної задачі 
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п-1 


К 
тах | з 3! арена 
тє| -1,1) кед Фум'обу- 
1 п-1 4 
Пе - Уси 4-» тіп 124 Х00. 
-1 каб Сомиоб ої 


Для них справедлива наступна теорема. 
Теорема 2. |9| При 4 2 со,1,2 поліномами найкращого наближення є відповідно 
поліноми Чебишова 1-го роду 
со5(п х агссо8 1) 
СА 


Т (І) - 


з 


поліноми Чебишова 2-го роду 


Та (1) - 


зіп((п- Т)агссозі ) 


пого 


а" 
(2п у а" 
Це означає, що при побудові поліноміальних інтерлінантів треба вибирати 


поліноми Лежандра 


ПО а С В 


прямі інтерлінації хх, і зт; узувіз Іп так, щоб числа хьізі,т; Увіз іп 
були коренями відповідних поліномів з найменшим відхиленням. Зокрема, при 4-1 


справедлива наступна теорема. 
ді Р2 


Теорема 3. (10) Нехай ре(рур,), Лу) є С"(/7), а|-1,1, р'з Ре 
дх |дх, 


вра 1802) | з є СР(72), рРя з 1) , ЕСД) - величина найкращого наближення функції / 


множиною В? за нормою |-||-| -| у 8 є В? - елемент найкращого наближення; 


а 
- со5(дл/(р, 1), йе1,рі, Хо, 7 СОБ(ЇЛСр» 1). із З 1,р, нулі поліномів 0, 


ХХ; 
їй 


Чебишова 2-го роду відповідно степеня р, та р»: 


Б р У ОД гу НИ 
зіп0 
ці і - базисні поліноми Лагранжа, ші, Са біо 
5 "(а У/ои о ій (хі) зн За р, із (2)- (1) 


й! і, зві 


РР. 
-2321 (ах Хо зір СА ро 0) 


й ЗТ, 1 
Для 709 існує єдиний, найближчий до 7(х) за нормою |||, елемент є є В? іцей 
елемент має вигляд (1), тобто є інтерлінантом, який інтерлінує /(х) на лініях 
Хо й зір, Каі1,2. 
Залишок наближення функції /(х) найкращим елементом має наступний вид. 


С, 20, 5) 


дртоРа 


то2-в'де зі в 2), бо) є 2 


рірз! 
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де (5,,6,) - деяка точка, що залежить від (х,,1ь)є/ 2, Тому найменше значення 


2 р 
величини || /(2)- 8 (2) досягається на тих 8 , для яких величина ||| || Ї (ху - Хі | є 
Ка зі ЗІ 
найменшою. Цю умову задовольняють поліноми, вузли яких є нулями поліномів 
Чебишова 2-го роду. 
2 Обчислення інтеграла від швидкоосцилюючих функцій двох змінних з викори- 
станням лагранжевої поліноміальної інтерлінації функцій. Для наближеного обчис- 


лення інтеграла 
ре ПР іанайноіьіі 
пропонується кубатурна формула ій 
уод- і банана зіпохахах», 
де ій 


ІЙ Оох)з Уои о ій са) зе р 027 ра і рі Арі і (ха). 


й! із, 1 1» З1 
В розгорнутому вигляді кубатурна формула має вид: 
Р Ха Ха ізн 
(7,9); Ї і й (х)зіпохах, Ї Ї одвіпохуд, чн 
4-1 Мій Хі 
р» Хорні Ман 
-У Ї Іра (0) біпохах, Ї Га, біпохах 
усі Хо Ма 
Р р» Мачі Хо із 
-У; 2. ху б? ГІ ій (х)зіпомах, Ї Ір, з (5 )віпо лох» . 
донос! Ма 2 


Теорема 4. Справедлива наступна оцінка похибки наближення інтеграла 1(/7, 0) 
кубатурною формулою 1 ( 7, с») 
М; 2 
Оо)- (ро) єс 
21 2(рі Її рь чі) 
Доведення. Знайдемо оцінку похибки наближення. Маємо 


1 о о)- Ї о) Пл зо )- (хо До з 


-1-1 


Р; Ро чані Зі 


С (м У, (25) 


РІТР2 
2 рірз! 


ПР кю, « 


дено зі Ха Ов 


Р; Ро чані Он 0, ад 


«мУУ | 


Ли Ха  Оі 


Р р» Ман іні 


2 р полин ра Ор а) дах, | 


з дїіьзі Хі 


на 


РІЯРа 
2 рірз 


ах, з 


0, , з Ж з 
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: ом р р- Ма зід(( ру -- 1) агосов х. ) ра х) ах чай набора зпоеано « 


РІ ні 1 2 
оле рірз 4-1, ЗІ й Тех Хі фе Х 
М У З Хі й ах Хо 4 ах, 


Зо 
РІтР2 | І 
2 руирз! 22 ) їз х 1 12- Р 
М Р; р» л Ха л 
2-1 3,2, агосовх, ---агосов х» 
Ма 2 


РІТР2 се 
2 рірз! дзії, сі 
М Р Р» 
регент 203 аГССО5 М 7 ЖГССО5 Ху; аГССО5 Хі, 7 агссо5 х,; 11Ї- 
РІтР2 1 2 2 
2 рірзі і -ьчі 


б М У» лі (і чі) лі - ті, 1) М 
2" руріцяьч РІЇ 0 рокі Дрокі | речі 
з М ВР п по Мрурот? | 
2" рірзїцль чі РІЗІре і 2Й"Р( р 1 рь НТ) 


Хоч 


і 


Теорема доведена. 
З Чисельний експеримент. Для функції 720,25) - с05(х, -25), у якої о у) - Її 


наведені в табл. | наближені значення інтегралів / ( й ,0) за кубатурною формулою 


(1,0) для різних значень р, - р, тафаелт, т. 1,2,.... 


Таблиця 1 - Обчислення І ( ї ,6) за кубатурною формулою Ї ( у ,6) 


о | рр; (17.9) 1(07о)- (09) ; 

Зл 6 4 -0.032615945938263 0.000000000292593 0.000000382471649 
6 -0.032615946229612 О.О00000000001244 0.000000003414925 
8 -0.032615946230859 ООООО00000000003 О.ОО000000001581 
9 -0.032615946230853 О.ОООО00000000003 ООООО00000000889 
10 -0.032615946230856 0 О.ОООО00000000045 

Зл 7 4 -0.032615945817734 0.000000000413122 0.000000027888558 
6 -0.0326159462291 О.ОО0000000001756 0.000000000249005 
8 -0.032615946230861 О.ОООО00000000004 О.ОО0000000001 153 
9 -0.032615946230852 ООООО00000000003 О.ОООО00000000065 
0 -0.032615946230856 0 ООООО00000000003 
4 -0.03261594623163 О.ООО000000000774 0.000000001770702 
6 -0.032615946230859 ООООО00000000003 О.ОО000000001581 
8 -0.032615946230856 0 О.ООО000000000073 
6 -0.018165043732214 О.ОО0000000000182 0.000000003414925 
7 -0.018165043732261 О.ОО0000000000135 О.О00000000013157 
8 -0.018165043732396 0 О.ООО000000000073 

рр 


Нехай є - . Наведемо точні значення інтегралів: 


2 (р, Ср, я 1)! 


1(7,3л) - -0.032615946230856 ; ЩО; Дл ) --0.018165043732396. 
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Висновки 


Побудована кубатурна формула наближеного обчислення інтеграла від швидкоос- 
цилюючих функцій двох змінних з використанням лагранжевої поліноміальної інтерлінації 


функцій на класі дійсних функцій, визначених на С - Г-1, її , і таких, що | ї (рр (х, у) «М. 


Отримана оцінка похибки кубатурної формули. Чисельний експеримент підтверджує те- 
оретичний результат. 
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Наближене обчислення подвійних інтегралів від швидкоосцилюючих функцій... 


О.М. Іуїууп, О.Р. МесПпиїуйег 

Арргохітате Саїсиатоп ої Доифіе Ппіеотаїз ої НіеПіу 
Озсійатогу Кипспопу улП Убе ої те Гаждтапое 
Роїупотіаї! ГпіегПпеапоп 


Модегп ргобіетея ої дїєта! віспа! ргосез5іпе пееад Ше з5о0Їийоп5 уліб пему Їогт58 ої даїа. 
І паеапя Шаг іпіогалайоп абоці Гипспоп 15 8еї ої уаїшез оп Ше 8еф ої Паїз ог 5еї ої Ппе5 ог 5еї 
ої Кпоїз. Тре Шеогу ої іпіегіїпеайопя апа іпіегНакайоп ої Типспопз5 15 Ше тобі еНесйує іп 
рів сазе. 

І Фе уогК, сибаштге Гогтиїа8 ої арргохітаї/е саїсціацоп ої аЧоцбіе іпіестаїя ої пієпІу 
о5сШаюогу Гипсйопбулті ц5е Фе І астапее роіупопиа! іпіегійпеацноп ої Гипсйопя а Ше сіа55 ої 
ФНегепбабіе Гипспоп5 аге іпуе5поаїед їп Фе агісіе. Пійогтлайоп абоші Гипсйоп 158 5еб ОЇ 
уаїше5 оп Ше орйтаї регрепдїсиціаг Ппез. Тре езштайоп5 ої епгог ої арргоасріпя, 0 Ше 
сибакиге Гогтиіаз аге ргезепіеа. 


Стаття надійшла до редакції 25.01.2012. 
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